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Abstract. Frequencies of 13 oscillation modes in the star 
XX Pyxidis (CD-24 7599) are accurately measured but 
for none of the modes the spherical harmonic degree (£) 
is known. We present results of an attempt to construct 
the model whose low-^ mode frequencies reproduce pos- 
sibly close the observations. Models are constrained by 
the mean photometric and spectroscopic data for the star. 
However, the strongest constraint on the effective temper- 
ture is from the requirement that the modes excited in the 
star fall into the range of the modes driven by the opacity 
mechanism. 

Our models are built with the standard stellar evo- 
lution code allowing no overshooting from the convective 
core. Effects of rotation are taken into account both in 
stellar evolution and in linear nonadiabatic oscillation cal- 
culations. Uniform rotation rate and conservation of the 
global angular momentum during evolution are assumed. 

We find several distinct mode identifications and as- 
sociated stellar models leading to frequency fits of simi- 
lar quality. Determination of the £ values for some of the 
modes could remove the ambiguity. None of the fits is sat- 
ifactory: the mean departures exceed the mean observa- 
tional frequency error by at least one order of magnitude. 
The fits could be improved by means of adjusting model 
parameters that were kept fixed. However, such effort will 
be meaningful only after improving accuracy in calcula- 
tion of the effects of rotation in oscillation frequencies. 

Key words: stars: oscillations - stars: evolution - stars: 
variables: S Scuti - stars: individual: CD-24 7599 = 
XX Pyx 



1. Introduction 

Long-time network observations of Delta Scuti stars re- 
sulted in establishing truly multimodal character of os- 
cillations in a few objects of this type. The four objects 
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with the largest number of modes detected are FG Virginis 
with 24 modes (Breger et al. 1997), 4 Canum Venaticorum 
with at least 17 modes (Breger 1997), XX Pyxidis with 
13 modes (Handler et al. 1997a) and BH Piscium with 13 
modes (Mantegazza et al. 1996). Though these numbers of 
modes in comparison with the Sun or certain oscillating 
white dwarfs may appear small, reproducing mode fre- 
quencies with model calculation represents a challenge. 

What makes our task more difficult than in the solar 
and white-dwarf cases, is the fact that the observed spec- 
tra do not reveal patterns enabling mode identification. 
Modes detected in S Scuti stars are of low radial order 
and their frequencies do not obey simple asymptotic re- 
lations. Further, these stars are rather rapid rotators and 
rotational splitting is not equidistant. Multicolor photom- 
etry may be used to determine the spherical harmonic 
degree, £. Spectroscopy allows, in addition, to determine 
the azimuthal order, to. In practice, however, reliable £ 
and to values are only seldom available. Thus, mode iden- 
tification in 8 Scuti stars cannot be done independently of 
calculations involving construction of equilibrium models 
and their oscillation properties. The aim is to construct a 
stellar model with frequencies of low-degree modes which 
fit the observed ones. In analogy to helioseismology we use 
the name seismic model. 

The ultimate purpose of seismic model construction is 
testing the ingredients of stellar evolution theory: its ba- 
sic assumptions and the input microscopic physics. Stellar 
evolution theory has reached an advanced level some 20 
to 30 years ago. However, important questions concerning 
convection and rotation remain unanswered. For instance, 
is the mixing- length theory sufficient for describing con- 
vective flux in subphotospheric layers? If so, what is the 
best choice for the mixing length parameter, a, depend- 
ing on stellar parameters? What is the extent of element 
penetration, d over , from the convective to the radiative re- 
gions? What is the law of the angular momentum evolu- 
tion? What is the role of rotationally induced instabilities 
in chemical element mixing? 
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In recent years considerable progress has been made in 
treatment of the equation of state and opacity for stellar 
interiors (see Christensen-Dalsgaard & Dappcn 1993 and 
references there, Rogers et al. 1996, Rogers & Iglesias 1994 
and references there, Iglesias & Rogers 1996, Seaton 1996). 
The uncertainties, however, are difficult to estimate. It 
is therefore important to use possibly diversified sets of 
observational data for testing. Observers provide us very 
accurate frequency measurements but a major effort is 
needed to connect these quantities to other stellar charac- 
teristics and to parameters used in model construction. 

In the next section we will outline methodology of con- 
structing seismic models, which, we believe, is applicable 
to all multimode 6 Scuti stars in the Main Sequence phase 
of evolution. We should stress that our work is not the 
first effort aimed at reproducing observed frequencies in 
objects of this type. Goupil et al. (1993) modelled GX Pe- 
gasi for which five frequencies have been measured. Breger 
et al. (1995) and Guzik and Bradley (1995) constructed 
models of FG Virginis, which approximately reproduced 
observations. Our target accuracy is much higher - we 
want to approach the level of the observational error. 

The rest of this paper concerns XX Pyxidis (CD- 
24 7599), which seemed to be a good choice for the 
first attempt in accurate reproducing observed spectrum. 
Amongst most multimode 5 Scuti stars it is the hottest 
and least evolved. In such stars the predicted spectrum of 
unstable modes is relatively simple. With the progress of 
stellar evolution the instability extends to mode of mixed 
pressure/gravity (p/g) character and to pure g-modes. 
Separations in frequency between modes of the same I > 
decrease, which complicates mode identifcation. In post 
Main Sequence stars the density of predicted spectra is 
so high that the task of mode identification will remain 
impossible unless we discover the clue to mode selection. 
Furthermore, the uncertainty in the description of the sub- 
photospcric convection is less severe in this relatively hot 
star. 

After presenting data in Section 3, we discuss in Sec- 
tion 4 constraints on mean parameters of the star. The 
methodology of the simultaneous search for mode identi- 
fication and refinement of model parameters is described 
in Section 5. In the same section several alternative solu- 
tions are presented. In Section 6 we explain why these are 
unreliable and in the last Section we discuss prospects for 
construction of a still improved seismic model of the star. 

2. Construction of seismic stellar models 

2.1. Principles 

The number of modes detected in individual <5 Scuti stars 
is far too small to attempt to determine radial structure 
directly from measured frequencies. Everything we may 
hope to achieve is to use the observed frequencies to de- 
termine global parameters characterizing the star and cer- 



tain parameters of the theory. The equations of problem 
are obtained by equating observed and calculated frequen- 
cies, 

/j\obs = fj, ca ,i(£j,mj,nj,P s ,PT), (1) 

where j labels measured frequencies, £j,mj 7 rij are num- 
bers identifying the mode, Ps gives the set parameters 
characterizing the model, and Pt gives the set of free pa- 
rameters of the theory. 

A comment is needed about rij. Nonradial modes ex- 
cited in S Scuti stars could be p-modes, g-modes, and 
modes of mixed character. In the present application, con- 
tinuity of / as function of P for specified (£, m, n) is essen- 
tial. The avoided crossing effect causes that the continuity 
is guaranteed if we define the n values at ZAMS where the 
g- and p-mode spectra are separated in frequency Follow- 
ing Unno et al. (1989) we will denote with n < and n > 
modes which in ZAMS models are consecutive g- and p- 
modes, respectively. For I > 1 there is the fundamental 
mode which is denoted with n = 0. With this choice, in 
evolved models some of modes denoted with n > will 
have prcdominatly g mode properties and vice versa. 

The parameters given by Ps are those characterizing 
the evolutionary sequence like initial mass M, chemical 
composition, Xq and Zq, and the angular momentum. In 
addition, there is a single parameter identifying the model 
within the evolutionary sequence. The age is always a good 
choice. In the present application we use the initial equa- 
torial velocity Kot,o instead of the angular momentum and 
log T ff instead of the age. The latter choice is acceptable 
because XX Pyx is certainly in the expansion phase of the 
Main Sequence evolution and T e g is a monotonic function 
of the age. In any case there are five parameters in Ps- 

The quantities one may include in Pt are mixing 
length parameter a, overshooting distance d ovel , as well as 
parameters characterizing angular momentum evolution 
and mass loss. None of these will be included in the present 
work. Constraining these parameters is certainly the most 
important application of asteroseismology. However, as we 
will see later in this paper, we are not yet at this stage. 

2.2. Tools 

2.2.1. Evolution code 

Stellar evolution code we use was developed in its original 
version by B. Paczyhski, R. Sienkiewicz and M. Kozlowski. 
It is in fact a modern version of B. Paczyriski's code 
(Paczyhski 1969, 1970), which is now written in a modu- 
lar form enabling an easy implementation of microscopic 
physics data from various sources. In the present applica- 
tion we use most recent versions of OPAL opacity (Iglesias 
& Rogers 1996) and equation of state (Rogers et al. 1996) 
data. The nuclear reaction rates are the same as used by 
Bahcall and Pinsonneault (1995). 
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The effect of averaged centrifugal force is taken into ac- 
count in the equation for hydrostatic equilibrium. We as- 
sume uniform rotation and conservation of global angular 
momentum. These are the simplest assumptions which we 
are prepared to abandon if this is required for a successful 
seismic model construction. Similarily, we are prepared to 
introduce overshooting and/or rotationally induced mix- 
ing of the chemical elements outside of the convective core, 
but these effects are ignored in the present work. With 
these assumptions the input parameters for sequences are 
total mass, M, initial values for hydrogen abundance, Xq, 
metal abundance, Zq, and the equatorial velocity Vrot,o- 
Initial heavy element mixture is that of Grevesse & Noels 
(1993). The assumed parameters for convection are a = 1, 
clover = 0. The choice of a is unimportant in the present 
application. 

For selected values of T e g the code returns all needed 
parameters of the model such as logg c ff, V~ ro t or fl and the 
matrix of the coefficients needed for calculation of linear 
nonadiabatic oscillation properties. 

2.2.2. Nonadiabatic oscillation code 

The code we use is a modified version of the code de- 
veloped long ago by one of us (Dziembowski 1977). The 
important modification is taking into account the effect 
of the averaged centrifugal force. Modified equations for 
adiabatic oscillations are given by Soufi et al. (1997). The 
corresponding change in the full nonadiabatic set is trivial. 
The input to this code specifies the range of frequencies 
and the maximum value of I. The standard output are the 
following characteristics of the m = modes: complex di- 
mensionless frequencies a, modified growth rates, 77, and 
the complex ratio of the flux to displacement eigenfunc- 
tions at the surface, T . 

The frequency, /, which we compare with observations 
is given by the real part of a: 

/ = ( 2 ) 

For higher order p-modes beginning with, say, n = 5 
or 6, the nonadiabatic effects are important at the level 
of the accuracy of frequency measurements. The growth 
rate is given by — Q c (c). A more convenient measure of 
mode stability is 77 which varies from —1, if damping oc- 
curs everywhere in the stellar interior, to +1, if driving 
occurs everywhere. At neutral stability we have rj = 0. 
Values of rj arc important as constraints on modes and 
models. We will use them in the next sections. Also JF's 
have applications in asteroseismology, especially, in deter- 
mination of ^-values (see Cugier et al. 1994). However, we 
do not have neccesary observational data on XX Pyx to 
make use of this quantity. 



2.2.3. Calculation of the rotational splitting 

Rotation has an important effect on the structure of oscil- 
lation spectra. Even at modest equatorial velocities such 
as 50-100 km/s the effect of rotation can not be reduced 
to the linear splitting oc mQ. Our code calculating the 
rotational splitting is accurate up to f2 2 . It is a version 
of the code by Dziembowski and Goode (1992) modified 
by two of us (WAD, AAP) and M.-J. Goupil. It uses the 
adiabatic approximation. The same nonadiabatic correc- 
tion evaluated with the the code described in the previous 
subsection is added to all modes within the multiplet. 

3. Frequency data 

More than 350 hours of time-series photometric observa- 
tions of XX Pyx are available. These include two Whole 
Earth Telescope runs (Handler et al. 1996, 1997a) as well 
as follow-up measurements to study the star's amplitude 
and frequency variations (Handler et al. 1997b). 

Analysis of the two WET data sets allowed the ex- 
traction of 13 pulsation frequencies plus one 2f-harmonic. 
These frequencies are summarized in Table [j], and mostly 
taken from Handler et al. (1997a, hereafter HPO). More- 
over, from an examination of the follow up data, the val- 
ues of frequencies f±, $2 and could be refined. The 
three revised frequencies are also listed in Table |l|. Their 
error sizes differ from those listed in HPO, since these 
frequencies are now known without alias ambiguity. How- 
ever, these three frequencies are slightly variable (Handler 
et al. 1997b). Therefore the error bars have been modified 
to account for the latter effect. 



Table 1. Frequencies of the 13 pulsation modes and the 
2/i -harmonic of the 5 Scuti star XX Pyx unambiguously de- 
tected by HPO 





Frequency 


Ampl. (1992) 


Ampl. (1994) 




(cycles/day) 


(mmag) 


(mmag) 


fx 


38.1101 


± 


0.0004 


11.5 ± 0.2 


15.9 ± 0.2 


h 


36.0113 


± 


0.0010 


10.0 ± 0.2 


3.3 ± 0.2 


h 


33.4370 


± 


0.0002 


6.4 ± 0.2 


3.8 ± 0.2 


h 


31.3925 


± 


0.0025 


3.8 ± 0.2 


1.9 ± 0.2 


h 


28.9950 


± 


0.0025 


2.4 ± 0.2 


2.7 ± 0.2 


h 


27.0028 


± 


0.0025 


2.4 ± 0.2 


1.6 ± 0.2 


fr 


34.6657 


± 


0.0038 


2.2 ± 0.2 


2.0 ± 0.2 


/b=2/i 


76.2203 


± 


0.0038 


0.8 ± 0.2 


1.1 ± 0.2 


h 


27.2635 


± 


0.0138 


1.0 ± 0.2 


1.9 ± 0.2 


ho 


29.6196 


± 


0.0075 


1.5 ± 0.2 


1.3 ± 0.2 


hi 


33.6367 


± 


0.0125 


1.2 ± 0.2 


1.7 ± 0.2 


fl2 


28.6980 


± 


0.0138 


1.1 ± 0.2 


1.5 ± 0.2 


fl3 


31.2033 


± 


0.0163 


1.4 ± 0.2 


1.1 ± 0.2 


/l4 


31.9062 


± 


0.0100 


1.0 ± 0.2 


1.7 ± 0.2 



4. Constraints 
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4--1- Mean surface parameters from observations 

To estimate the position of XX Pyx in the H-R diagram, 
color photometry and high-resolution spectroscopy of the 
star have been carried out by HPO. Applying calibra- 
tions to the photometric data, they infer that XX Pyx is a 
main-sequence S Scuti star with T e g — 8300 ± 200 K, log 
g = 4.25±0.15. From their spectroscopy HPO determined 
the star's projected rotational velocity with v sin i — 52±2 
km/s and that the object has approximately solar metal 
abundance: [M/H] = 0.0 ± 0.2. 

4-2. Structures in the power spectra 

An important constraint for identifying the different pul- 
sation frequencies of XX Pyx with the corresponding 
quantum numbers was found by HPO. They searched for 
characteristic spacings within the 13 pulsation frequen- 
cies with a Fourier technique: they assumed unit ampli- 
tude for all the individual pulsations they detected and 
Fourier analysed the resulting signal. This is analogous to 
a spectral window: the power spectrum of such a signal 
has maxima at the sampling frequency and its harmonics 
(see Kurtz 1983). With this method it is possible to find 
regular frequency spacings in an objective way. 

There are two physical reasons why regular frequency 
spacings can be present: in case of slow (rigid) rotation 
the members of different multiplets are approximately 
equally spaced in frequency. Since the measured v sin i 
of XX Pyx is 52 km/s (HPO), rotational splitting could, 
under favourable circumstances, be detected. The second 
possibility is that consecutive modes of the same I are ex- 
cited. Even for low-overtone pulsations such modes are 
approximately equally spaced in frequency. 

To illustrate this second possibility, we used our basic 
grid of the main-sequence models and their oscillation fre- 
quencies (see subsection 5.2): M — 1.75 — 2.05 M with 
a step of 0.05 M©, logT cff = 3.905 - 3.925 with a step 
of 0.05, Vrot.o = 50 — 125 km/s with a step of 25 km/s. 
For each of these models, we compared the asymptotic fre- 
quency spacing with the mean spacing between frequen- 
cies of consecutive modes of a given £ value. We restrict 
to the observed frequency range of XX Pyx. The results 
are shown in Fig. [jj. (Remember that the asymptotic fre- 
quency spacing for high overtone p-modes is determined 
by the following formula, see Unno et al. 1989, Eq. (16.36): 
Ai> — 0.5/ (J Q =£), where c is the sound speed, — i.e. it is 
proportional to the square root of the star's mean density.) 

As can easily be seen from Fig. ??, the assumption that 
regular frequency spacings are present even at relatively 
low overtones is justified. The frequency spacings of the 
different ts are very similar. Consequently, HPO carried 
out the Fourier analysis described above and discovered 
a regular frequency spacing within the pulsation modes 
of XX Pyx, with a probability of less than 3% that this 
detection was caused by chance. Because of the value of 
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Fig. 1. Asymptotic frequency spacing of high-order p-modes 
of our main-sequence models compared to the mean spacing of 
1 = 0,1,2 modes in the frequency range excited in XX Pyx. 
Only axisymmetric (m = 0) modes are shown. Rightmost 
points correspond to stellar models of lowest mass and highest 
effective temperature. The higher mass and/or the lower effec- 
tive temperature the higher overtones of p-modes are appeared 
in a fixed frequency range, therefore the mean frequency spac- 
ing is approaching the asymptotic value. The scatter of the 
datapoints, most pronounced for I — modes, reflects the dif- 
ferent rotati on rate s of our models. The discontinuities of I = 2 
modes at ^ (p/po) < 0.46 are due to intruding g-modes, while 
more regular level shifts of £ = 1 and £ = 2 modes are caused 
by different overtones (and different number of these overtones) 
appearing in the frequency range of interest. 

this spacing (« 26 /iHz w 2.25 c/d), HPO suggested that it 
cannot reflect rotational splitting (since the second-order 
terms would destroy any frequency symmetry at such high 
rotation rates). 

From a comparison of the Fourier analyses of the ob- 
served and of several model frequency spectra, HPO con- 
cluded that this spacing is caused by the presence of 
modes corresponding to consecutive overtones of t = 1 
and I = 2 modes. HPO determined the spacing of con- 
secutive overtones of the same £ with 54.0 ± 2.3 /iHz 
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(= 4.67 ± 0.20 c/d). This frequency spacing is a measure 
of the sound crossing time through the object, and thus it 
can be used to determine the mean density of the star. For 
a model sequence with solar metallicity, HPO determined 
p = 0.246 ± 0.020 Pq. 

HPO performed several tests to verify the validity of 
this method and the correct interpretation of its result 
(see their paper for more details) and found their Fourier 
technique to be reliable, since the frequency spectrum of 
XX Pyx is very well suited for such an analysis: only a 
few overtones are excited, but a large number of possible 
pulsation modes are photometrically observed. 

Actually, one can use this Fourier method to test a 
common assumption in mode identification attempts from 
photometric (and radial velocity) data: with increasing de- 
gree £ of the pulsations, geometrical cancellation (Dziem- 
bowski 1977, Goupil et al. 1996) decreases the photomet- 
ric amplitude of the modes. Therefore, it is usually as- 
sumed that only modes with £ < 3 can be detected using 
photometric data. The photometric amplitude of an £ = 3 
mode with the same intrinsic amplitude as an £ = 1 mode 
is only about a factor of 12 smaller, which is approximately 
the range of photometric amplitudes of the different pul- 
sation modes detected for XX Pyx. However, when con- 
sidering rotational splitting, seven modes of £ = 3 can be 
present, but only three £ = 1 modes. Therefore it can not 
be ruled out that a non-negligible number of £ = 3 modes 
is excited to observable amplitude. 

Consequently, we applied HPO's Fourier technique to 
model frequencies, incorporating modes with £ — 3 and/or 
£ = 4. In cases where many of these high £ modes were 
present, no significant frequency spacing could be found, 
since their rotationally split patterns "masked" the regular 
spacing of the I — — 2 modes. Only with a few of the high 
£ modes the typical frequency spacing could be revealed. 
This suggests that most of the pulsation modes of XX Pyx 
are indeed of £ = — 2. 

4-3. Implication from the stability survey 

XX Pyx is a hot S Scuti star located near the blue 
edge of the instability strip. In such an object the opacity 
mechanism driving is confined to modes in a narrow range 
of frequencies and the driving effect is marginal. This is 
seen in Figure |[ Values of r\ never exceed 0.05. The range 
of / where r\ > never spans more than three radial orders 
for radial modes. In models considered the lowest values of 
n is 4 and the highest is 8. The behaviour of T](f) in several 
models is shown in Fig. 0. The selected models have the 
surface parameters in the range consistent with the data 
on XX Pyx. The independence on £ and occurence of a 
single maximum are typical features for the opacity driven 
modes. 

Values of rj do not have the same reliability as those 
of /. Our treatment of the optically thin layers and the 
convective flux is crude. Therefore we will regard identi- 
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Fig. 2. Normalized growth rates, n, plotted against frequency, 
/, in selected models spanning the allowed ranges of M and 
logTgff. The values of M (in solar units) and logT^ are given 
in each panel. Models were calculated with Xo = 0.7, Zq — 0.02 
and Vrot = 0. The symbol ps is plotted near the corresponding 
radial overtone. The thick horizontal line shows the range of 
frequencies measured in XX Pyx. 

fications with modes of small negative 77's as admissible. 
However, we can exclude identifications involving radial 
modes with n < 4 and n > 8 as well we can exclude all 
models with, approximately, logT ff > 3.925 and the low 
mass models with logT c ff > 3.92 (see Fig. |||). 

5. A search for best parameters and mode identi- 
fication 

We have no observational information about the £ and m 
values for the modes detected in XX Pyx. The assignment 
of these two quantum numbers as well as n may only be 
done together with model parameter determination. The 
basis is Eq. (1) applied to the 13 measured frequencies. 
We only assume that all modes are of low degree. Most 
of the results we will present here were obtained with the 
assumption £ < 2 for all modes. We carried out also cal- 
culation allowing one of the mode to have t — 3. Some of 
the results for such a case will be discussed at the end of 
this section. 

At this stage we allow no free theoretical parameters 
Pt and we fix Xq and Zq at the values 0.7 and 0.02, re- 
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spectively. Thus, we allow only three free model parame- 
tres: M, V I0 t,o, and T c g. For each model the modes are 
assigned to minimize 

X ~j '^2/i—li.fohs.i /calc,i) ) (3) 

where j — 13 or 12. The first value was used in the case 
when we assumed I < 2 for all modes. In the alternative 
version we search for the minimum of x 2 relaxing the fit 
for one of measured frequency. In the automatic fitting 
routine we made sure that each of the model frequencies 
was assigned at most to one of the stellar frequencies. 

5.1. Uncertainties in stellar and model frequencies 

In our search for best models and mode identification we 
use unweighted observational data. The reason is that the 
errors in the frequency measurements are smaller than 
the uncertainities in the theoretical values. In Table [j] we 
may see that the errors range from 4 x 10~ 4 for f\ to 
163 x 10~ 4 for /13. Thus the relative accuracy of frequency 
is, at worst, ~ 5 x 10 , which is by one order of magnitude 
less than the theoretical uncertainty. 

The main problem on the side of theory is in the 
treatment of rotation. We found that at equatorial veloci- 
ties ~ 100 km/s the second order perturbation treatment 
yields relative accuracy in frequencies ~5x 10 -3 . 

5. 2. Tabulation of model frequencies 

With the codes described in subsection 2.2 we prepared 
tables with frequencies for modes with I = 0, 1, 2 cover- 
ing the range of the stellar pulsation frequencies. At this 
stage we neglected the effects of near-degeneracy which 
we discuss in Sect. 6, because it complicates calculations 
and it is not essential at this stage. The range of logT c ff 
values was [3.905 — 3.925] with a step of 0.05. The lower 
limit is somewhat less than the lower limit allowed by 
the photometry (subsection 4.1) while the upper limit fol- 
lows from the stability consideration (subsection 4.3). The 
range of equatorial velocities was [50 — 125] km/s with a 
step of 25 km/s. Here the lower limit follows from the 
usini measurement (subsection 4.1). The upper limit was 
adopted to avoid large errors in treatment of rotation. The 
adopted mass range was [1.75 — 2.05]M Q with a step of 
0.05M Q . The implied range of the mean density is sig- 
nificantly larger than that inferred in subsection 4.2. We 
consider models in the wider range of mean densities be- 
cause there is a nonzero probability that the HPO spacing 
is caused by chance. 

5.3. Properties of X 2 {M, V Iot , , T e $ ) 

We calculated values of \ 2 according to Eq.(3) with j = 13 
for nearly 40000 models interpolating model parameters 
and frequencies from the three-dimensional basic grid de- 
scribed above. At this stage we confined the upper limit of 



2 - 




Fig. 3. Values of x 2 evaluated for frequencies in ~ 7000 models 
in the ranges M = 1.75 - 2.05 M Q , Kot,o = 50 - 110 km/s, 
logTcff = 3.905 — 3.925. The patterns are the same (but less 
dense) as in the case of ~ 40000 models used to choose best 
identifications. 

the equatorial rotational velocity still more, to the value 
of 110 km/s. In Fig. || we plot the values of \ 2 against 
mean density, which is the most important parameter 
determining p-mode spectra. We see several dips in \ 2 
which allows us to isolate potential mode identification. 
One of the dips occurs close to p = 0.246, which was the 
preferred value from the observational frequency spacing. 
The present analysis shows that there are few alternative 
"good" values. We stress that a good fit may occur even 
if there is no preferred spacing. The coincidence however 
speaks in favor of identfications corresponding to this par- 
ticular dip. The dips become shallower with the increasing 
p which merely reflects the fact that in less evolved models 
there are fewer modes in the specified frequency range. 

At this stage we have to consider the identifications 
associated with the dips located near the following eight 
values of p: 0.17, 0.19, 0.22, 0.25, 0.30, 0.36, 0.41 and 0.49. 
It should be noted, however, that all the minimum values 
of x 2 are still far larger than the measurement errors. 

Next important parameter is the equatorial velocity of 
rotation which determines the m-dependence. The plot of 
X 2 against V To t,0, which we do not reproduce here, shows 
two shallow minima around the values 60 and 100 km/s. 

5.4- Possible models and identifications 

There is no doubt that if all (or almost all) 13 modes 
detected in XX Pyx are of I < 2 degree then the mode 
identification and model must correspond to one of the 
dips in x 2 ■ At each dip we chose a number of identifica- 
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tions with nearly the same lowest values of \ 2 - All such 
identifications are listed in Table ||. They are grouped ac- 
cording to the dips. For each specified identification the 
three model parameters - M, V lot .o, logoff - are deter- 
mined through the \ 2 minimalization. These and other 
model characteristics are listed in Table pi 

Not all models listed in Table || are consistent with the 
requirement that the identified modes are pulsationally 
unstable. Thus far we made a limited use of the constraint 
following from the stability considerations. 

In fact, the stability argument may be used to elimi- 
nate some of the models listed in Tables || and || as can- 
didates for the seismic model of XX Pyx. A comparison 
of the model positions with the blue edges in the H-R 
diagram, shown in Fig. [|, facilitates the elimination. We 
see that all the models with M < I.QMq and hotter than 
some logT e ff = 3.91 lie rather far from the blue edge for 
modes with / < 27 c/d, which is the lowest frequency 
observed in XX Pyx. At higher masses the upper limit 
of log T c ff should be moved to 3.914. In this way we may 
eliminate models D2B, D2C, D3B, D3C and D8A. 

Perhaps model DlA may be eliminated with the oppo- 
site argument because it lies far to the red of the 27 c/d 
blue edge. This means that there are unstable modes in 
the model with lower frequencies that are not seen in the 
star. This is, however, a much weaker argument because 
it is not an uncommon situation that unstable modes are 
not excited with detectable amplitudes. 

In Fig. H the frequencies determined for XX Pyx are 
compared with frequencies of low degree modes in selected 
models. These are not all acceptable but they cover the 
whole range. The common pattern of the model frequency 
spectra is the large departure from the equidistant pattern 
of the rotational splitting. This is particularly well seen in 
the case of the £ — 1 triplets. One may see that with the in- 
creasing radial order the triplet becomes more and more 
asymmetric, so that the prograde mode m = — 1 almost 
overlap the centroid mode m = at the high frequency 
end. There is a near equidistant separation between the 
consecutive centroid modes at each £ in most of the cases. 
The exception is the low frequency part of the £ = 2 spec- 
trum in all, but the lowest mass models. Here we see the 
effect of avoided crossing between two multiplets. Modes 
in this range have mixed (p- and g-) character. Their fre- 
quencies are sensitive to the structure of the deep interior 
and therefore these modes are of special interest for seis- 
mic sounding. 

5.5. Fitting 12 modes 

Naturally, if we demand fitting for only 12 of the measured 
frequencies we may attain significantly lower values of \ 2 ■ 
To estimate the effect we looked for the minimum of x 2 
considering all possible choices of 12 modes for each model. 
As expected, the general patterns of the \ 2 dependence on 
pIPq were found very much like those in Fig. ^. In partic- 




Fig. 4. Position (symbols) of the best models in the H-R dia- 
gram. The asterisks show the models used in Fig. |^. The evo- 
lutionary tracks for models with indicated mass and V TO t = 
are shown with dashed lines. The line denoted with BE shows 
the absolute blue edge; there are no unstable modes in mod- 
els to the left of this line. The lines denoted with p n are the 
blue edges for the corresponding radial overtones and the line 
denoted with 27 c/d is the blue edge for modes with the cor- 
responding frequency. 

ular, the dips occured at the same locations. For the same 
set of models the absolute minimum of \ 2 (0.00254) was 
reached at dip 3 (p/p© — 0.216). The minimum value was 
nearly three times smaller than the absolute minimum in 
the previous 13-mode case. The implied mean frequency 
mismatch of 0.05 c/d is close to the estimated uncertainty 
in calculated frequencies. It is, thus, clear that improve- 
ment on the side of theory is needed before we will be able 
to produce a credible seismic model for this star. In the 
next section we discuss two most needed improvements. 

6. Problems 

One important effect in the treatment of rotational fre- 
quency perturbation which we did not include so far in 
our calculation was the near-degeneracy of certain modes 
which may be coupled by rotation. The effect is discussed 
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Table 2. Possible identifications of the 13 modes detected in XX Pyx. Horizontal lines separate groups of the identifications 
associated with the consecutive dips of \ 2 seen in Fig. ^. Note that within the group most of the identifications are identical. 
Models corresponding to each identification are determined by minimalization of \ 2 ■ 
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h 


/l2 
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/io 


/is 


h 


/l4 


h 


/n 


fr 


h 


h 




27.003 


27.264 


28.698 


28.995 


29.620 


31.203 


31.392 


31.906 


33.435 


33.637 


34.666 


36.012 


38.110 


Model 


1 n m 


1 n m 


1 n m 


1 n m 


1 n m 


1 n m 


1 n m 


1 n m 


1 n m 


1 n m 


1 n m 


1 n m 


1 n m 


D1A 


2 3-1 


2 3-2 


15 1 


2 4-1 


1 5-1 


2 5 


2 5-1 


2 5-2 


1 6 


1 6-1 


2 6 1 


2 6-2 


2 7 2 


DIB 


2 3-1 


2 3-2 


1 5 1 


2 4-2 


1 5-1 


2 5 


6 


2 5-2 


1 6 


1 6-1 


2 6 1 


2 6-2 


2 7 2 


D1C 


2 3-1 


2 4 2 


1 5 1 


2 4-2 


1 5-1 


2 5 


2 5-1 


2 5-2 


1 6 


1 6-1 


2 6 1 


2 6-2 


2 7 2 


DID 


2 3-2 


2 4 2 


1 5 1 


2 4-2 


1 5-1 


2 5 


2 5-1 


2 5-2 


1 6 


1 6-1 


2 6 1 


2 6-2 


2 7 2 


DIE 


2 3-1 


2 4 2 


1 5 1 


2 4-2 


1 5-1 


2 5 


6 


2 5-2 


1 6 


1 6-1 


2 6 1 


2 6-2 


2 7 2 


D1F 


2 3-2 


2 4 2 


1 5 1 


2 4-2 


1 5-1 


2 5 


6 


2 5-2 


1 6 


1 6-1 


2 6 1 


2 6-2 


2 7 2 



D2A 
D2B 
D2C 


1 4-1 
1 4-1 
1 4-1 


2 3-1 
2 3-2 
2 3-2 


2 4 
2 4 
2 4 


5 
5 
2 4-1 


1 5 1 

2 4-2 
2 4-2 


1 5 
1 5 
1 5 


1 5-1 
1 5-1 
1 5-1 


2 5 1 
2 5 1 
2 5 2 


2 5-1 
2 5-1 
6 


1 6 1 

1 6 1 

2 5-1 


2 6 2 
2 6 2 
1 6 1 


2 6 1 
2 6 1 
2 6 2 


2 6-2 
2 6-2 
2 6-2 


D3A 
D3B 
D3C 


2 3 
2 3-1 
2 3-1 


1 4 1 
14 1 
14 1 


2 4 2 
14 
2 4 2 


1 4-1 
1 4-1 
1 4-1 


2 4 1 
2 4 1 
2 4 1 


5 
5 
5 


2 4-1 
2 4-1 
2 4-1 


1 5 1 
1 5 1 
1 5 1 


1 5 
1 5 
1 5 


1 5-1 
1 5-1 
1 5-1 


2 5 1 
2 5 1 
2 5 1 


2 5-1 
2 5-1 
2 5-1 


1 6 
1 6 
1 6 


D4A 
D4B 


2 3 2 
2 3 2 


2 3 1 
2 3 1 


4 
2 3-1 


2 3-2 
2 3-2 


14 1 
14 1 


1 4-1 
1 4-1 


2 4 2 
2 4 2 


2 4 1 
2 4 1 


5 
5 


2 4-2 
2 4-2 


1 5 1 
1 5 1 


1 5-1 
1 5-1 


6 
6 


D5A 


1 3 1 


2 2-2 


1 3-1 


2 3 2 


2 3 1 


2 3-1 


4 


2 3-2 


14 


2 4 2 


2 4 1 


2 4 


1 5 1 


D6A 


2 2 2 


2 2 1 


2 2-1 


2 2-2 


1 3 1 


1 3 


1 3-1 


2 3 2 


2 3 


2 3-1 


2 3-2 


1 4 1 


2 4 2 


D7A 
D7B 


2 1-2 
2 1-2 


1 2-1 
1 2-1 


2 2 2 
2 2 2 


2 2 1 
2 2 1 


2 2 
2 2 


3 
3 


2 2-2 
2 2-2 


1 3 1 
1 3 1 


1 3 
1 3-1 


1 3-1 

2 3 2 


2 3 1 
2 3 1 


2 3 
2 3 


14 1 
14 1 


D8A 


2 10 


2 


2 1-2 


1 2 


1 2-1 


2 2 2 


2 2 1 


2 2 


2 2-2 


3 


1 3 1 


1 3 


2 3 1 



Table 3. Parameters of models determined by the \ 2 minimalization. Nm is the total number of modes with I < 2 within 
the frequency range of modes detected in XX Pyx. The number decreases with mean density which explains increase of the \ 2 
minima. Note that the difference in mean radial mode degree between consecutive groups is close to 0.5. Models used in Fig. g 
are marked with asterisks in the last column. 



Model 


M/Mq 


Kot(ZAMS) 




logT cff 


logL 


p/Pe 


x 2 


N m 


E(n + //2)/13 




D1A 


1.9869 


57.50 


51.92 


3.90533 


1.2906 


.16692 


.006876 


30 


5.92 




DIB 


2.0181 


61.82 


55.92 


3.91101 


1.3175 


.16717 


.006217 


29 


5.92 




D1C 


2.0160 


60.18 


54.43 


3.91063 


1.3158 


.16710 


.006224 


29 


6.00 




DID 


2.0266 


59.38 


53.69 


3.91247 


1.3249 


.16708 


.007508 


29 


6.00 




DIE 


2.0165 


61.55 


55.69 


3.91072 


1.3161 


.16715 


.006108 


29 


6.00 


* 


D1F 


2.0260 


60.59 


54.79 


3.91235 


1.3243 


.16711 


.007646 


29 


6.00 




D2A 


1.9953 


108.01 


99.30 


3.91241 


1.2876 


.18680 


.008070 


27 


5.73 


* 


D2B 


2.0477 


110.35 


101.56 


3.92169 


1.3320 


.18681 


.007450 


26 


5.62 




D2C 


2.0333 


66.67 


61.20 


3.92056 


1.3233 


.18842 


.007387 


26 


5.62 




D3A 


1.9122 


95.00 


88.65 


3.90567 


1.2064 


.21608 


.007856 


24 


5.23 


* 


D3B 


1.9986 


97.64 


91.16 


3.92154 


1.2819 


.21614 


.007181 


23 


5.23 




D3C 


1.9972 


95.40 


89.09 


3.92132 


1.2809 


.21605 


.007530 


23 


5.23 




D4A 


1.8722 


66.67 


63.00 


3.90611 


1.1635 


.24711 


.019895 


20 


4.77 




D4B 


1.8782 


67.27 


63.60 


3.90723 


1.1688 


.24732 


.020005 


20 


4.77 


* 


D5A 


1.8367 


81.32 


78.62 


3.90907 


1.1149 


.29872 


.036899 


19 


4.23 


* 


D6A 


1.8527 


60.28 


59.07 


3.92234 


1.1148 


.36149 


.045757 


18 


3.69 




D7A 


1.7556 


75.00 


74.63 


3.90944 


1.0074 


.41550 


.050614 


17 


3.27 




D7B 


1.7562 


74.13 


73.72 


3.90933 


1.0086 


.41368 


.048300 


17 


3.31 




D8A 


1.7500 


61.11 


61.17 


3.91500 


0.9819 


.48689 


.061828 


18 


2.85 
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M = 2 016, <p/p B > = 0,167, log T e „ = 3,911, log L = 1.316, V rol = 56 km/s 
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Fig. 5. Comparison of the XX Pyx frequencies (vertical lines) 
with model frequencies (symbols) for £ = (top row), 1 = 1 
(middle row), 1 = 2 (bottom row). The encircled symbols de- 
note the modes identified with those excited in the star. Small 
open circles denote the stable modes (always very close to in- 
stability 77 ~ 0. Positions of the selected models in the H-R 
diagram are denoted with asterisks in Fig. In Tables |2| and ^ 
these models (from top to bottom) are denoted as DIE, D2A, 
D3A, D4B, and D5A. 



in details by Soufi et al. (1997). Here we will present se- 
lected results for our model D3A. 

Uniform (more generally, spherical) rotation couples 
modes if their azimuthal numbers (m) are the same and 
their degrees (£) are the same or differ by 2. If the fre- 
quency difference between the two (or more) coupled 
modes is of the order of the rotation frequency one has 
to use the version of the peturbation formalism appropri- 
ate for the case of degeneracy i.e. to consider as a zeroth 
order basis a linear combination of the nearly degener- 
ate modes. The relative contributions of the components 
are determined from the second order (in Q) perturbation 
equations. 

In application to XX Pyx a systematic near-degene- 
racy occurs between the £ = and 2 modes and be- 
tween the £ = 1 and 3 modes. In the latter case, for each 
multiplet we have three coupled modes corresponding to 
m = —1, 0, 1. These near-degeneracies follow from rather 
high values of n and, hence, approximate validity of the 
p-mode asymptotics. We also considered coupling involv- 
ing three modes, e.g. £ = 0, 2, 4 or £ = 0, 2, 2, in the case 
of avoided crossing. In none of the cases considered the 
inclusion of the third mode was essential. 

Model D3A and its oscillation frequencies were ob- 
tained by means of interpolation described in the previous 
section. In order to evaluate the effect of the degeneracy we 
had to recalculate the model and its frequencies. A com- 
parison of the two upper panels of Fig. || clearly shows 
that the grid of models was not sufficiently dense. The 
differences between interpolated and calculated frequen- 
cies are easily seen, especially for the £ = 2 modes near 
the avoided crossing. This inadequacy of our grid is not 
essential at this stage but must be kept in mind in more 
advanced efforts. In the top panel one may see how the fit 
is improved if the 1 = 2 identification for = 34.7 c/d is 
replaced with the £ — 3. The value of \ 2 is then lowered 
by some 30 %. No fit improvement is achieved by allowing 
the I — 3 identification for another poorly fit frequency 
/io = 29.6 c/d. One may see, however, in the mid panel 
that the situation changes if instead of the interpolated 
frequencies one uses the calculated ones. 

A comparison of the mid and bottom panels shows the 
effect of the coupling between nearly degenerate modes. 
The effect is best seen in the case of I — 2 and modes. 
The frequency distance between the modes increases when 
the coupling is taken into account. In Table |^ we give the 
values of the frequency shifts (Sf) caused by the mode 
coupling as well as the amplitudes of the spherical har- 
monic components (A\ ow is corresponding to i = or 1 
and Ahigh is corresponding to £ — 2 or 3). One may see 
that the mutual contamination of the £ = and 2 com- 
ponents is quite strong. This is bad news for prospects of 
mode discrimination by means of two-color photometry 
because such contaminated modes may appear in the am- 
plitude ratio - phase difference diagrams (Watson 1988, 
Garrido et al. 1990) in the £ = 0, 1 or 2 domains depend- 
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ing on the inclination of the rotation axis. Fortunately, 
modes with m^O are not affected. The mutual contam- 
ination is also significant for the £ = 1 and 3 pairs, fn 
all the cases the nominal £ = 3 modes will be most likely 
observable through their £ = 1 contaminations. 



Table 4. Effects of coupling between nearly degenerate modes 



26 28 



30 32 34 

Frequency [c/d] 



I 

38 



40 



Fig. 6. Effects of inadequate grid density and near degener- 
acy shown for one of the best models (D3A). The upper panel 
shows interpolated frequencies which are the same (except that 
£ = 3 modes are added) as in the mid panel in Fig. ^. The mid 
panel here shows frequencies of the same modes for the recal- 
culated model. Frequencies of the modes coupled by rotation 
are connected with dashed lines. The resulting frequencies of 
these modes are shown in the lower panel. 



7. Discussion 

Clearly, we have not succeeded in constructing the seismic 
model of XX Pyx. Models regarded as plausible, such as 
those in Fig. |[ do not reproduce the frequencies within 
the observational errors. In fact, we are quite far from 
the 1CP 3 — 10~ 2 c/d error range seen in Table @. The de- 
partures from the fit for some of the identifications are in 
fact easily visible in Fig. ^. The problem proved more diffi- 
cult than we have anticipated. We believe that we learned 
something in the process and that this knowledge should 
be shared with other groups undertaking similar efforts. 

Undoubtedly, an improvement of the fits could be ac- 
complished by allowing adjustment of the chemical com- 
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-0.385 
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27.239 


0.112 


0.935 


-0.356 
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3 


1 


26.720 


-0.112 


0.422 


0.906 


fs 


1 


4 


-1 


28.932 


0.096 


0.951 


-0.308 




3 


3 


-1 


28.293 


-0.096 


0.374 


0.927 


/l3 





5 





31.153 


0.226 


0.881 


-0.474 
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4 





30.678 


-0.226 


0.513 


0.858 


/l4 
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31.923 


0.153 


0.915 


-0.403 
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4 
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31.443 


-0.153 


0.481 


0.877 


fs 
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5 





33.374 


0.122 


0.962 


-0.275 
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4 





32.302 


-0.122 


0.337 


0.942 


hi 


1 


5 


-1 


33.613 


0.137 


0.931 


-0.365 




3 


4 


-1 


33.038 


-0.137 


0.441 


0.897 


h 


1 


6 





38.064 


0.174 


0.944 


-0.329 




3 


5 





36.988 


-0.174 


0.370 


0.929 



position parameters. The overshooting distance and pa- 
rameters in the radial dependence of rotation rate should 
also be regarded as adjustable quantities. These degrees of 
freedom affect primarily the positions of the mixed modes 
relative to pure p-modes. In this application the mixed 
modes exist for £ = 2 at the low frequency end. Such 
modes extend to higher frequencies for higher ts. How- 
ever, the freedom allows also a fine tuning of the distances 
between pure p-modes. One could hope that there is just 
one mode identification and one corresponding model of 
the star for which the frequency fit within the measure- 
ment errors is possible. 

The reason why we cannot yet proceed this way is that 
our treatment of rotation is not adequate. We used here a 
perturbation theory which is accurate up to some 0.05 c/d, 
much worse than the measurement accuracy. We are not 
yet sure that the recently developed cubic theory (Soufi 
et al. 1997) would yield a sufficient accuracy. 

Observational determination of the £ values for some of 
the excited modes could significantly change the situation. 
Imagine, for instance, that we find that the photometric 
data place the f\ mode into the I = domain of an ampli- 
tude ratio - phase difference diagram (see Watson 1988, 
Garrido et al. 1990). We may see in Table [l] that this mode 
has far the highest amplitude and certainly it will be the 
first for which we will have the required data. As we have 
discussed in Section 6, the possibility that this is not a 
£ — mode but rather a £ = 2, m = mode must be con- 
sidered. Fortunately, we do not see such an identification 
for /i in Table |^. Thus, we are left with only two models, 
D4A and D4B, with very similar parameters and the same 
identification for almost all modes excited in XX Pyx. We 
could discriminate between the two if we knew £ for fi2 
which is the only mode with different identification. If, 
however, the f\ position in the diagram corresponds to an 
£ = 1 domain, we are less lucky. The possible models are 
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not only six models with I = 1 for /i but also models D4A 
and D4B with I = because radial modes contaminated 
with a i = 2 component may appear in the I — 1 domain. 
The ambiguity implies a wide range of admissible model 
parameters. 
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